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ABSTRACT 

We present the first measurement of the mass function of galaxy clusters based directly on cluster masses derived 
\^ ' from observations of weak gravitational lensing. To investigate the degree of sample incompleteness resulting from the 
. X-ray based selection of the target clusters, we use a sample of 50 clusters with weak lensing mass measurements to 
■ empirically determine the relation between lensing mass and X-ray luminosity and the scatter about this relation. We 
' use a complete, volume-limited sub-sample of 35 X-ray luminous clusters of galaxies at 0.15 < z < 0.3 to constrain 
^jTj: the abundance of very massive (Af > l0^^h~^MQ) clusters. From this, we constrain o'8(rim/0.3)°'^^ — 0.671q (68% 
^ confidence limits), agreeing well with constraints from the 3-year WMAP CMB measurements and estimates of cluster 
, abundances based on X-ray observations, but somewhat lower than constraints from "cosmic shear" weak lensing 
^ • measurements in random fields. 

CO ' Subject headings: Cosmology: observations — dark matter — gravitational lensing — large-scale 

CN i structure of the Universe — galaxies: clusters 
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1. INTRODUCTION 

It has long been recognized that the abundance of mas- 
sive clusters of galaxies is a sensitive probe of density fluc- 
tuations in the universe. The cluster mass function can 
be used to derive interesting constraints on parameters 
that govern structure growth, in particular the normal- 
ization of the primordial spectrum of density fluctuations 
and the average mass density of the universe Qm (e.g., 
Evrard 1989; Frcnk et al. 1990; Henry & Arnaud 1991; 
Lilje 1992; White, Efstathiou, & Frenk 1993; Cen & Os- 
triker 1994; Eke, Cole, & Frenk 1996; Henry 1997; Ouk- 
bir & Blanchard 1997; Bahcall & Fan 1998; Pen 1998; 
Viana & Liddle 1999; Henry 2000). 

If the matter-energy density in the universe is domi- 
nated by a dark energy component, structure growth is 
hindered when this component becomes dynamically im- 
portant. The resulting effect on the evolution of the clus- 
ter mass function makes it an interesting probe of dark 
energy (e.g., Wang & Steinhardt 1998; Haiman, Mohr, 
& Holder 2001; Levine, Schulz, & White 2002; Munshi, 
Porciani, & Wang 2004; Henry 2004; Majumdar & Mohr 
2004). Being a sensitive probe of the growth factor, clus- 
ter abundances might be combined with probes of the 
expansion history such as the magnitude-redshift rela- 
tion for type la supernovae to constrain photon-axion 
oscillations (e.g., Csaki, Kaloper, & Terning 2005) or 
even modifications to gravity arising in brane-world cos- 
mologies (Lue, Scoccimarro, & Starkman 2004). Finally, 
cluster abundances can provide limits on the level of 
non-Gaussianity in the primordial density perturbations 
(Scoccimarro, Sefusatti, & Zaldarriaga 2004; Mathis, 
Diego, & Silk 2004). 

^ Based on observations made with the Nordic Optical Tele- 
scope, operated on the island of La Palma jointly by Denmark, 
Finland, Iceland, Norway, and Sweden, in the Spanish Observato- 
rio del Roque de los Muchachos of the Instituto de Astrofi'sica de 
Canarias 

^ Visiting observer, University of Hawaii 2.24m Telescope at 
Mauna Kea Observatory, Institute for Astronomy, University of 
Hawaii 



Many current studies use the X-ray temperature func- 
tion (XTF) or the X-ray luminosity function (XLF) of 
galaxy clusters as proxies for the mass function (e.g., 
Pierpaoh et al. 2001; Ikebe et al. 2002; Reiprich & 
Bohringer 2002; Seljak 2002; Allen et al. 2003; Pierpaoh 
et al. 2003; Viana et al. 2003; Henry 2004; Stanek et al. 
2006). Assuming that the X-ray emitting intra-cluster 
medium (ICM) is in a state of hydrostatic equilibrium, 
the observed X-ray temperatures Tx and luminosities Lx 
can be simply related to cluster masses. However, this 
assumption is not always valid, as some clusters have 
complex temperature structures due to ongoing merging 
events (e.g., Henry & Briel 1995; Henriksen & Marke- 
vitch 1996; Pratt, Bohringer, & Finoguenov 2005). A 
significant spatial offset between the peaks of the X-ray 
emission from the ICM and the corresponding concentra- 
tions of (collision-less) dark matter is observed in some 
of these systems (Clowe, Gonzalez, & Markevitch 2004; 
Jee et al. 2005). Simulations indicate that major merg- 
ers will temporarily boost ICM temperatures and lumi- 
nosities significantly above their equilibrium values and 
thereby significantly alter the observed XTF and XLF 
(e.g., Randall, Sarazin, & Ricker 2002). This could po- 
tentially bias constraints on cosmological parameters, if 
the resulting scatter around the XTF and XLF is not pre- 
cisely determined for clusters in the high-mass tail of the 
distribution (i.e., if the box sizes of the simulations are 
too small to adequately sample the most massive clus- 
ters). 

These problems can be remedied by obtaining indepen- 
dent values for the cluster masses through measurements 
of gravitational lensing, which do not rely on assump- 
tions about the dynamical state of the clusters. A set of 
lensing-derived cluster masses can then be used to em- 
pirically calibrate the normalization and scatter in the 
M - Lx and M - Tx relations (e.g., Allen et al. 2003; 
Smith et al. 2005; Pedersen & Dahle 2006). StiU, the 
complicated gas processes in merging clusters may plau- 
sibly lead to a significantly asymmetric dispersion about 
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the mean relations which can not be adequately probed 
by the limited current data sets of clusters with both 
good X-ray and lensing data. 

A different strategy is to derive the cluster mass func- 
tion directly from lensing masses. The main advantage of 
this approach is that it should in principle avoid the prob- 
lem of relating any observables of the baryonic cluster 
component to actual cluster mass (in practice, however, 
such issues are still relevant whenever the cluster sam- 
ple is selected based on baryonic tracers of mass, but at 
a less problematic level) . This approach requires lensing 
masses for a large sample of clusters with well- understood 
selection criteria. Samples from cluster surveys that ef- 
ficiently pick out the extremely rare, most massive clus- 
ters in a large volume are most suitable, since the high- 
mass end of the mass function is most sensitive to cos- 
mological parameters. Also, given the current level of 
statistical and systematic uncertainties in lensing mea- 
surements, very massive systems are needed to produce 
a useful lensing signal from individual clusters. How- 
ever, the lensing measurements are subject to projection 
effects that provide an additional source of error, and 
need to be calibrated to reasonable precision, e.g. using 
N-body simulations (e.g., Clowe, de Lucia, & King 2004). 

Here, we present the first attempt to constrain the clus- 
ter mass function by using weak gravitational lens mea- 
surements directly. Our results are only dependent on 
observables of the baryonic component of the clusters in 
the sense that the clusters were selected based on their 
X-ray luminosity. We describe our sample selection in 
§ 12 and in § O we describe in detail how we use our 
mass measurements to derive the mass function. Partic- 
ular emphasis is put on the treatment of statistical and 
systematic errors and how these affect our results. 

Throughout this paper we will assume a spatially-flat 
cosmology with fi,„ -I- fl\ — 1. Numbers are generally 
given for the "concordance model" case (ftm, ^a) = 
(0.3, 0.7). The Hubble parameter is generally given by 
Hq ~ lOO/i kms~^Mpc~^ (however, in the discussion of 
the M — ix relation we use i/o = 5G/i5o kms^^Mpc^^ 
for easy comparison to results in the literature). 

2. CLUSTER SAMPLE SELECTION 

Cluster surveys based on photometric and spectro- 
scopic followup of the ROSAT All-Sky Survey (RASS; 
Triimper 1993) using optical telescopes are currently the 
most efficient strategy for constructing representative 
samples of very massive clusters. These surveys cover 
large volumes, but are not restricted to optically selected 
objects such as the clusters in the catalogs of Abell (1958) 
and Abell, Corwin, & Olowin (1989), which arc seriously 
incomplete in the redshift interval relevant for this work 
(Briel & Henry 1993). 

The cluster sample used here was selected from a larger 
sample of X-ray luminous clusters of galaxies which were 
imaged using monolithic 2048^ CCD detectors and/or 
the 8192^ UH8K mosaic CCD camera. We obtained 
V- and /-band imaging data for the clusters using the 
2.56m Nordic Optical Telescope on La Palma, Canary 
Islands, Spain and the 2.24m University of Hawaii Tele- 
scope at Mauna Kea Observatory, Hawaii, USA. Expo- 
sure times were typically 1.5h in each passband with the 
2048^ detectors and 3.5h in each passband with the less 
sensitive UH8K camera. The seeing was in the range 



O'.'e < FWHM < I'.'l in both pass-bands, with a mean 
seeing of 0'.'81 and 0'.'89 in the 7-band and F-band, re- 
spectively. For the weak lensing analysis, background 
galaxies were selected based on signal to noise ratio 
rather than magnitude, the limits corresponding to the 
magnitude ranges 21 < m/ < 24.5 and 22 < my ^ 25.5 
for point sources. The typical number density of lensed 
sources was ~ 25 galaxies per square arcminute, and 
we found a "figure of merit" value (defined by Kaiser 
2000) of Y.Q Id^ - 1-5 X lO^deg-^. Further details 
about the observations, data reduction and weak lens- 
ing analysis are given elsewhere (Dahle et al. 2002; H. 
Dahle 2006, in preparation). The reduction procedures 
for mosaic CCD data followed the techniques described 
by Kaiser et al. (1999). The clusters used here repre- 
sent a volume-limited sample of X-ray luminous clusters 
taken from the RASS-based, X-ray flux hmited ROSAT 
Brightest Cluster Sample (BCS) of Ebehng et al. (1998) 
and its low-flux extension (eBCS; Ebeling et al. 2000). 
We will hereafter refer to them collectively as the (e)BCS 
sample. We chose a lower cutoff in X-ray luminosity 
of -^x,o.i-2.4kcV = 6 X lO^'' erg s^^ (limit given for a 
concordance model universe with h = 0.7). As dis- 
cussed in § 13.31 and illustrated in Figure ^ this cutoff 
approximately corresponds to a lower cutoff in mass of 
Afisoc 5 X lO^/i-^Afg. Following White (2002), we 
denote as Migpc the mass contained inside the cluster ra- 
dius risoc within which the average density is 180 times 
the critical density of the universe Pcr:{z) at the redshift 
of the cluster. 

2.1. Sample Completeness 

The estimated overall completeness of the BCS sample 
within its X-ray flux limit is 90%, and the completeness 
of the eBCS sample is estimated to be 75% (Ebeling 
et al. 1998; 2000). The selected clusters lie within the 
redshift range 0.15 < z < 0.303. There are 35 (e)BCS 
clusters that fall within the given limits (and also do 
not have their X-ray signal signiflcantly contaminated by 
point sources), and we have made weak lensing measure- 
ments of all of these (see Table Another cluster with 
a high lensing mass, CI 1821-^643 (Wold et al. 2002), 
is located within the survey volume, but it was omitted 
from the (e)BCS samples because of the presence of a 
strong AGN point source in the cluster. Bohringer et 
al. (2001) find the occurrence of such AGN point source 
contamination in a large X-ray selected cluster sample at 
similar redshifts to be less than 5%, so we do not expect 
to miss more than 1-2 additional clusters in our sample 
from this effect. Furthermore, such cases are accounted 
for in the BCS and eBCS completeness estimates quoted 
above, and we therefore chose not to include a weak lens- 
ing mass measurement for CI 1821-1-643 in our sample. 

Ebeling et al. (1998; 2000) identified the majority of 
their (e)BCS clusters from the data base produced by 
running the Standard Software Analysis System (SASS) 
on the RASS data (Voges 1992). The SASS data have 
a total exposure time of 360s with the ROSAT Position 
Sensitive Proportional Counter (PSPC), this exposure 
time being constant across the sky. Photon event table 
(PET) files were extracted in a 2° x 2° region centered 
on each of the SASS-based cluster candidates. In some 
cases (particularly near the ecliptic poles) the PET fields 
(which contain all the RASS photons collected in the re- 
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gion) are significantly deeper than the SASS data, which 
led to serendipitous discoveries of additional clusters us- 
ing a Voronoi tessellation and percolation algorithm (for 
full details, see Ebehng et al. 1998). 

The BCS and eBCS completeness values quoted above 
are average values; for both samples there is a small mi- 
nority of clusters that are serendipitous detections, and 
have much lower completeness values associated with 
them. This is the case for two clusters in our sample 
(A1576 and ZW5247). For each cluster, we hst the es- 
timated completenesses, calculated from the complete- 
ness associated with the relevant (e)BCS cluster detec- 
tion method, in Tabled 

3. CLUSTER MASS FUNCTION 

A number of steps, detailed below, are required to con- 
vert weak gravitational lensing measurements for a set 
of clusters into an observed mass function which can be 
compared to predictions from theory or N-body simula- 
tions. Firstly, the total survey volume is calculated for a 
range of cosmologies, and lensing estimates of the cluster 
mass within a spherical volume are derived. The effects 
of sample incompleteness should be taken into account, 
arising both from incompletenesses in our original clus- 
ter catalog (discussed in the preceding section) and the 
mass-dependent incompleteness caused by the observed 
scatter around the Migoc — -^x relation, coupled with 
applying a sharp lower cutoff in Lx- Another impor- 
tant correction comes from uncertainties in the cluster 
mass estimates, which will tend to artificially increase 
the number of clusters in the high-mass tail of the ob- 
served mass function. 

3.1. Survey volume 

The (e)BCS survey area consists of the Northern celes- 
tial hemisphere (6 > 0°), minus a "zone of avoidance" at 
Galactic latitudes |6| < 20°, where high Galactic extinc- 
tion may seriously decrease the depth of foUowup op- 
tical extragalactic observations. The survey volume is 
bounded in redshift by 0.15 < z < 0.303. For the case 
of the concordance cosmological model, this represents a 
volume of 8.0 x 10*(/i"iMpc)3. 

3.2. Cluster mass estimates 

For our weak lensing analysis, we used the shear esti- 
mator of Kaiser (2000). This method was tested, along 
with several other shear estimators, by Heymans et al. 
(2006), using simulated weak lensing data. The shear 
estimator of Kaiser (2000) is more mathematically rig- 
orous than the currently most widely used shear esti- 
mator (Kaiser, Squires, & Broadhurst 1995), but it dis- 
plays a significant non-linear response to shear, unlike 
most other shear estimators. If we correct our measure- 
ments using a second order polynomial in the difference 
between measured and true shear, as found by Heymans 
et al. (2006), we find that most cluster masses stay within 
+/ — 15% of the mass calculated from uncorrected shear 
values. Furthermore, the change in average cluster mass 
is < 2%, i.e., there is very little systematic shift in mass. 
Given the modest size of the correction, and the fact 
that it would (for a few clusters) require extrapolations 
outside the range of shear values over which the shear 
estimator has been tested, we chose not to apply this 
correction. 



The galaxy shape distortions caused by gravitational 
lensing provides a measurement of the reduced tangen- 
tial shear, gx = 7r/(l — k), where 7t is the tangen- 
tial component of the shear and k is the convergence. 
The average value of gx is measured in a set of non- 
overlapping annuli with mean radius r, centered on the 
mass center of the cluster. We fit a NFW-type pro- 
file (Navarro, Frenk, & White 1997) to the observed 
shear profile gxir) of each cluster. In this model, the 

density is p{r) oc [{r/rs){l + r/rs)'^'\ , and it is spec- 
ified by a scale radius r^ and a concentration param- 
eter Cvir = r^ir/rs, where rvir is the virial radius of 
the cluster. N-body simulations of clusters in a ACDM 
universe predict the dependency of Cvir on cluster mass 
and redshift (Bullock et al. 2001; Dolag et al. 2004). 
We have previously shown that the NFW model pro- 
vides a good fit to the average shear profile of a sam- 
ple of 6 galaxy clusters at z ~ 0.3 (Dahle, Hannestad, 
& Sommer-Larsen 2003), with a concentration parame- 
ter consistent with predictions. Studies based on X-ray 
observations reach similar conclusions (Pointecouteau, 
Arnaud, & Pratt 2005; Vikhhnin et al. 2006; Voigt & 
Fabian 2006). As shown in Tabled the median cluster 
mass of our sample is Migoc — 8 x 10^** M0, for which 
Bullock et al. (2001) predict a median halo concentra- 
tion of Cvir = r^iv/rs = l.lArisoc/rs = 5.8/(1 -|- z) in 
a concordance model universe. Given the modest range 
in masses for the clusters studied here, the weak mass- 
dependence (Cvir oc M~°-^^) predicted by these authors 
is negligible, and a fixed value for Cvir is adopted at a 
given redshift. For derivations of the lensing properties 
of the NFW model, we refer the reader to the papers of 
Bartelmann (1996) and Wright & Brainerd (2000). 

The NFW model fit gives estimates of rigoc and Afigoc. 
The shear measurements used for the fit were made at 
cluster-centric radii 50" < r < 180" for clusters which 
were observed with 2048^ CCD cameras and 150" < r < 
550" for the clusters which were observed with the UH8K 
mosaic CCD camera. For our clusters, we derive rigoc 
values typically corresponding to an angular extent in the 
range 450" < rigoc < 750". Hence, extrapolations were 
in many cases required to derive the masses within rigoc- 
As we assume a fixed value for Cvir at a given redshift, any 
intrinsic scatter in this parameter will introduce an extra 
uncertainty in the cluster mass estimates. The level of 
dispersion predicted by Bullock et al. (2001) and Dolag 
et al. (2004) (a la scatter of A(logCvir) ^ 0.18) gives 
a contribution to the overall error budget for the clus- 
ter masses which depends on the ratio of the maximum 
radius of the shear measurements r^t, to r'lsoc which is 
listed in Table ^ For the clusters observed with 2048^ 
CCD cameras (rfit/''i80c — 0.3), the predicted scatter 
in Cvir results in an additional mass uncertainty of 18%, 
which was assumed to be unrelated to other measure- 
ment uncertainties. Similarly, for the clusters observed 
with the UH8K camera {rfn/Tisoc — 1-0), an uncertainty 
of 10% was estimated and added in quadrature to other 
contributions to the mass measurement uncertainty (see 
below). 

Our photometric data in two pass-bands did not al- 
low reliable discrimination between cluster galaxies and 
lensed background galaxies (except for the small frac- 
tion of galaxies that had observed V — I colors redder 
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TABLE 1 

Weak lensing masses and X-ray luminosities 
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22 
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4.66 


3.37 


0.29 


0.90 


A586 


0.17 


13.0 


26.80 


8.89 


0.17 


0.90 


A1914 


0.17 


21.6 


7.20 


4.41 


0.26 


0.90 


A665 


0.18 


19.3 


8.27 


4.38 


0.26 


0.90 


A115 


0.20 


6.9 


6.01 


4.14 


0.31 


0.90 


A520 


0.20 


17.4 


13.03 


4.48 


0.75 


0.90 


A963 


0.21 


12.6 


6.65 


5.81 


0.95 


0.90 


A1423 


0.21 


12.1 


11.80 


5.94 


0.26 


0.90 


A773 


0.22 


16.0 


12.07 


5.38 


0.27 


0.90 


A2261 


22 
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1 
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A267 


0.23 


16.9 


11.80 


4.36 


0.86 


0.90 


A1682 


0.23 


13.8 


4.38 


2.92 


0.39 


0.90 


A1763 


0.23 


18.4 


8.34 


4.02 


0.31 


0.90 


A2111 


0.23 


13.4 


4.88 


2.85 


0.38 


0.90 


A2219 
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0.23 


26.4 


11.99 


5.28 


0.28 


0.90 
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0.23 


12.4 
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0.47 


0.15 


RX 791 9Q 


0.23 


22.9 


5.97 


4.09 


0.36 


0.90 
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5.08 


0.31 
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A1835 


0.25 


48.2 


8.39 


4.43 


0.34 


0.90 


A68 


0.26 


18.7 


17.08 


9.01 


0.27 


0.90 


M^^l 4'^'^-^99 


0.26 


16.6 


10.16 


4.71 


0.33 


0.90 


Zw 5768 


0.27 


14.7 


6.04 


4.65 


0.40 


0.78 


A697 


0.28 


20.8 


20.40 


7.34 


0.28 


0.90 


A1758N 


0.28 


14.9 


21.02 


7.36 


0.27 


0.90 


A2631 


0.28 


16.8 


4.88 


3.33 


0.44 


0.78 


A611 


0.29 


17.4 


5.21 


3.47 


0.29 


0.78 


RX J0437.H-0043 


0.29 


15.7 


3.74 


2.65 


0.49 


0.78 


Zw 3146 


0.29 


34.0 


8.94 


4.23 


0.38 


0.90 


Zw 7215 


0.29 


14.4 


8.38 


4.67 


0.38 


0.78 


A781 


0.30 


22.2 


10.00 


4.55 


0.37 


0.78 


A1576 


0.30 


14.1 


14.49 


4.66 


1.00 


0.13 


A2552 


0.30 


19.6 


3.70 


2.91 


0.52 


0.78 



Note. — The luminosity and lensing mass values are given for a flat cosmology with Q-m = 0.3, Ha = 0.7. The tabulated mass 
uncertainties include an assumed 26% from projection effects, added in quadrature to the lensing mass estimate. The tabulated errors are 
calculated by symmetrising the error bars. 

— See Ebeling et al. (1998; 2000) for references to redshift measurements. 
^ — The X-ray luminosity values are taJsen from Ebeling et al. (1998; 2000). 



than early- type cluster galaxies at the cluster redshift). 
The faint galaxy catalogs which were used to measure 
the gravitational shear would thus be significantly con- 
taminated by cluster galaxies that will dilute the lensing 
signal by an amount given by their local sky density. 
To correct for this contamination, a radially dependent 
correction factor was applied to the shear. The magni- 
tude of the correction was estimated based on the radial 
dependence of the average faint galaxy density in two dif- 
ferent sets of clusters observed with the UH8K camera, 
one at z ~ 0.3 and the other at z ~ 0.2, assuming that 
the contamination is negligible at the edge of the UH8K 
fields. This procedure led to corrections to the cluster 
masses at the 20-30% level (see Pcderscn & Dahle 2006 
for details). Clustcr-to-cluster variations in their rich- 
ness in faint dwarf galaxies will lead to an additional 
scatter when a fixed level of contamination is assumed 
for all clusters. Based on our data, we estimate a typical 
scatter in cluster dwarf richness of 50%, leading to an 
additional uncertainty in the cluster mass measurements 
of 12% for data obtained with with the UH8K camera 
(which include lensing measurements at large radii where 



the contamination is less severe) and 20% for the clus- 
ters observed with the 2048^ CCD cameras. These were 
added in quadrature to the mass measurement uncertain- 
tics. 

The redshift distribution of the background galaxies 
needs to be known to convert the observed lensing sig- 
nal into a cluster mass estimate. The background galaxy 
redshifts were estimated from spectroscopic and photo- 
metric rcdshifts in the Hubble Deep Field (for details, 
see Dahle et al. 2002). The average value of the ratio 
between the lens-source and observer-source angular di- 
ameter distances, /? = Dis/Ds, was approximated by a 
relation on the form = Az^^ + Bzc\ + C, where Zc\ is 
the cluster redshift. The constants were calculated for a 
range of spatially-flat cosmologies; for the concordance 
model, we obtained A = 1.37, B = -2.00, and C = 1.01. 

Gravitational lensing effects measure the projected 
mass distribution in the cluster, including mass outside 
''180c- This may introduce a bias and a dispersion in 
the derived lensing estimates of Migoc. Based on N- 
body simulations of clusters in a ACDM model, Met- 
zler, White & Loken (2001) considered a lensing mass 
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estimator based on the projected mass within r2oo6 and 
'"500& (the radii within which the mean density contrast 
is i5 = 200 and 5 = 500, respectively). For their estima- 
tor, they found that lensing masses would on average be 
biased by a factor (Afions/^true) = 1-33, with a mass dis- 
persion of 0.26 around the mean. However, these authors 
remarked that this positive bias could largely be removed 
by assuming a realistic model for the radial mass distri- 
bution out to large radii. A more recent study (Clowe, 
De Lucia, & King 2004) indicated that there is on aver- 
age no significant bias in the lensing-derived masses, pro- 
vided that these are determined by fitting the observed 
tangential shear to an NFW model, using shear measure- 
ments all the way out to Risoc- Still, a significant scatter 
will be added to the derived lensing masses, and we added 
such a scatter in quadrature to the measurement uncer- 
tainties of our lensing mass estimates. We performed 
our calculations using three different values of the scat- 
ter (0.26, 0.13, and 0), but found only a weak effect on 
the end result, as the shear measurement uncertainties 
dominate the error budget for most of the clusters. All 
subsequent results in this paper are based on adopting 
the highest value (26%) of this scatter. 

3.3. Mass-X ray luminosity relation 

The observed scatter around the Misoc — Lx relation 
causes a mass-dependent sample incompleteness which is 
most severe at the lower end of the range of measured 
cluster masses. Our sample is defined by a fixed cutoff in 
X-ray luminosity which corresponds to ix,o.i-2.4keV — 
1.2 X W^^h^Q erg s~^ for the concordance model. Any 
intrinsic scatter around the mean Migoc - Lx relation 
and uncertainties in the X-ray luminosity measurements 
tabulated by Ebeling et al. (1998; 2000) will contribute 
to soften the corresponding cutoff in mass. To assess the 
completeness of our sample as a function of Migoc, we 
made an empirical determination of the Migoc - Lx re- 
lation and its scatter. This was done using a data set of 
50 clusters (Dahle et al. 2002; H. Dahle 2006, in prepa- 
ration) which, in addition to the 35 clusters used here 
for the determination of the mass function, include clus- 
ters of similarly high X-ray luminosities, selected from 
four cluster samples (Ebehng et al. 1996, 1998, 2000; 
Bohringer et al. 2000), all based on RASS data. We 
adopted the RASS-based X-ray luminosities listed by 
these authors (with preference for the (e)BCS values 
whenever a cluster is included in more than one sample) , 
after scaling them to our chosen cosmological model. 

For many of these clusters, more precise Lx values ex- 
ist in the literature, e.g., from pointed observations with 
ROSAT, XMM-Newton or Chandra. However, our cur- 
rent aim is not to determine the Migoc — Lx with the 
highest possible accuracy, but rather to quantify the in- 
completeness effects caused by selecting clusters based on 
their X-ray luminosities in the RASS cluster surveys. Lx 
only serves as a proxy for mass for our cluster sample se- 
lection, and when determining (below) the scatter about 
the mean Migoc — Lx-relation we do not distinguish the 
contribution to the scatter caused by uncertainties in the 
Lx-nieasurements from the intrinsic scatter about the re- 
lation, as they will both act to soften the low-mass cutoff 
in a similar way. 

The lensing masses and X-ray luminosities of the clus- 



ters are plotted in Figure ^ The large apparent spread 
in this figure, compared to similar diagrams from X-ray 
data, is mostly a result of the narrow range in Lx, and 
hence, in cluster mass. 
We assumed a Migoc — Lx relation on the form 



E{z)Misoc = Mo 



Lx 



E{z) 



(1) 



where E{z) ^ [I + z)^{l + zQ^ + 17a/(1 + z)2 - Q^) 
is an evolution parameter that describes the redshift- 
dependence predicted from simulations (e.g., Bryan & 
Norman 1998; Mathiesen & Evrard 2001). Using the 
BCES(X2|Xi) estimator of Akritas & Bershady (1996), 
wc fitted the data in Figure ^ to a model on the form 



log 



E{z)M: 



180c 



a log 



Lx 



i;(z)io44/ir^'"--=-i 



log 



o 
Mo 



This estimator takes into account the errors in both 
axes, and allows for possible intrinsic scatter in the 
Misoc — Lx relation. Our best-fit values are a = 1.04 ± 
0.46 and log (A/o/Zi^Mq) = 14.0 ± 0.6. As can be ex- 
pected from Figure ^ the slope a is rather poorly con- 
strained, since the clusters in our sample only span mod- 
est ranges in Migoc and Lx- Figure ^ also seems to in- 
dicate a flattening of the Misoc — Lx slope at very high 
X-ray luminosities (ix,o.i-2.4kcV ^ 2 x 10'*^/i^o^ erg s~^), 
similar to the change of the slope of the M — Tx relation 
at high temperatures noted by Cypriano et al. (2004) . In 
both cases, this is most likely caused by a scattering of 
out-of-equilibrium clusters into the extreme tails of the 
XTF and XLF, as discussed in §□ 

We also tried fixing the slope parameter at the value 
a = 0.75 predicted from a self-similar collapse model for 
galaxy clusters (Kaiser 1986). Our fit to the relation 
in Eq. ^then gave a normalization log(Mo//i^Q^AfQ) = 
14.34 ± 0.07. This is consistent with the relation sim- 
ilarly derived from a sample of 17 clusters by Allen et 
al. (2003). These authors used a combination of X-ray 
based mass measurements derived from Chandra data for 
a set of 10 clusters they had classified as relaxed systems, 
and weak lensing-based mass measurements for a par- 
tially overlapping sample of 10 clusters, including some 
unrelaxed systems. Since Lx is our only cluster selection 
criterion, we have not here attempted to classify clusters 
according to their dynamical state. 

When accounting for the additional scatter about this 
relation caused by the mass measurement uncertainties, 
a scatter of (Tm — 0.44 is found. By comparison, the data 
presented by Reiprich & Bohringer (2002) indicate an in- 
trinsic scatter of ctm = 0.39 about the M — Lx relation, 
using X-ray based mass estimates under the assumptions 
of isothermality and hydrostatic equilibrium of the intra- 
cluster gas. More recently, Smith et al. (2005) found an 
intrinsic scatter ctm = 0.41 about their M — Lx rela- 
tion, based on gravitational lensing measurements of the 
masses of 10 clusters within a cluster-centric radius of 
250ft.~^ kpc. The slightly larger scatter we measure could 
be a result of the larger uncertainties of Lx values de- 
termined from RASS data, compared to deeper pointed 
X-ray observations of clusters. 



(2) 
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Fig. 1. — Weak lensing mass estimates of 50 clusters plotted against their X-ray luminosities. The dashed line indicates the best fit to 
the mass-luminosity relation, with an arbitrary slope a. Also shown is the best-fit relation when assuming the theoretically expected value 
a = 0.75 (solid line). The dotted lines represent the Icr uncertainty about this relation. 



For a given cluster of true mass Migoc, the probability 
of being included in our sample is 

PsmAM - M{L^in))dM, (3) 

-oo 

where Psml is the probability distribution of Afigoc for 
a given X-ray luminosity and M(Lmin) is the mass cor- 
responding (via eq. ^ to the cutoff in X-ray luminosity 
^min = 1-2 X lO^^h'^Q erg s~^. Here, Psml is assumed to 
be a log-normal function (analogous to the scatter about 
the X-ray luminosity-temperature relation; see Novicki, 
Sornig, & Henry 2002), for which we derive a standard 
deviation in AlogAf of ciogM = 0.178, based on the 
data plotted in Figure ^ We note that the exact form 
of Psml should be determined empirically in the future, 
using even larger data sets of X-ray luminous clusters 
with more accurately determined lensing masses. 

To make our subsequent analysis self-consistent, the 
appropriate value of M{Lrain) was determined for the 
range of spatially-flat cosmologies probed in this study; 
its value at the mean cluster redshift (z = 0.23) for the 
concordance model is M(Linin) = 5.2 x 

3.4. Measurement uncertainties 

A significant uncertainty is associated with any re- 
alistic cluster mass estimator. In the case of gravita- 
tional Icnsing-based mass measurements, uncertainties 
arise both from the measurement uncertainties of the 
observable itself and from the intrinsic scatter in the 
mass-observable relation caused by projection effects. As 
discussed in § 13.21 the former is dominant for most of 
the clusters in our data sample, but there is also a non- 
negligible contribution from the latter. Given these un- 
certainties, extra care must be taken to remove possible 
biases produced by the exponential falloff of the mass 
function at high masses. If the uncertainties in Migoc 
were the same for all the clusters, they could simply be 
incorporated by convolving the theoretical mass function 
with a probability distribution P(6Misoc), with a width 



given by the uncertainties (e.g., Evrard 1989; Metzler et 
al. 2001). However, as shown in Tabled the uncertain- 
ties vary significantly from cluster to cluster. Hence, the 
theoretical mass function is convolved with the uncer- 
tainty of each cluster, and the average of the resulting 
ensemble of mass functions is calculated and used in our 
further analysis. 

Given a theoretical mass function n.th(.^^i80c), and 
given the total measurement uncertainties in the weak 
lensing estimates of Afigoc, the observed mass function 
can be predicted as: 

nohs{Misoc) = 

Q dSMisocP{SMi80c)nth{Misoc + SMi8nc)P{Mi8ac)^ , 

where P(6Misqc) is a log-normal probability distribution 
of width given by the measurement uncertainty in A/igoc, 
and the average is taken over the ensemble of observed 
uncertainties ctmisoc- 

Finally, we note that Misoc, ctmisoc and M(imin) 
are all cosmology-dependent quantities, and to make our 
analysis self-consistent, these need to be re-calculated for 
each cosmology. 

3.5. Fit to predicted mass functions 

A number of studies have predicted the mass func- 
tion of dark matter halos using an analytic approach or 
via numerical simulations. The formalism of Press & 
Schechter (1974) uses linear theory for the evolution of 
self- gravitating mass density fluctuations and assumes a 
Gaussian distribution of inhomogeneities to calculate the 
fraction of the volume which has reached a given over- 
density for a given smoothing scale, assuming that all ob- 
jects at a spherical over-density above a threshold value 
Sc = 1.69 undergo gravitational collapse. The mass func- 
tion derived via Press-Schechter theory provides overall 
a reasonable match to mass functions derived from N- 
body simulations of CDM particles, over a wide range of 
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M18O 

Fig. 2. — Observed cumulative mass function. The apparent flattening of the slope at < 8x W^'^h"^ Mq is caused by severe incompleteness 
of the sample at lower masses. This p lot is shown for illustrational purposes only; the actual procedure used for fitting the data to theoretical 
mass functions is described in S I3.5I 



masses (e.g., White et al. 1993). However, the mass func- 
tions derived from simulations still show significant de- 
partures from the Press-Schechter form, yielding a larger 
number of clusters at the high-mass end (e.g., Jenkins 
et al. 2001). Sheth, Mo & Tormen (2001) showed that a 
modification of Press-Schechter theory to allow for ellip- 
soidal, rather than spherical, collapse gave a significantly 
improved fit to the simulations, similar to the form sug- 
gested by Sheth & Tormen (1999). Recently, Warren 
et al. (2006) used a new, larger set of simulations and 
found that the mass functions introduced by Jenkins et 
al. (2001) and Sheth & Tormen (1999) deviate by > 30% 
from their simulations at the highest masses. They there- 
fore provide a new formula for the mass function shape. 

Based on these results, we chose to fit our cluster data 
to two mass function forms; one which is consistent with 
the most recent simulations (Warren et al. 2006), and one 
which is more directly based on analytical work (Sheth 
& Tormen 1999). Jenkins et al. (2001) and White (2002) 
showed that the shape of the mass function is univer- 
sal across the range of currently interesting cosmologies, 
given that the halo mass (Migob) is always defined within 
a volume of average over-density S = 180 with respect to 
the mean (rather than critical) matter density of the uni- 
verse. In the case of the simulations. Warren et al. (2006) 
remark that the commonly used "friends-of-friends" halo 
definition method tends to select halos of a different av- 
erage over-density, which depends on the number of par- 
ticles in the simulated halo. For the most massive halos, 
the halo mass is M2sab, i-e., corresponding to an aver- 
age halo over-density A = 280. We convert the mass 
function predictions of Warren et al. (2006) and Sheth & 
Tormen (1999) to Migoc (from M28O6 and Migob, respec- 
tively), using the fitting function in Appendix C of Hu 
& Kravtsov (2003), based on the procedure outlined by 
White (2001). We assume a structure parameter fixed at 
r = 0.21. 

For easy comparison between theory and observations, 
we binned the theoretical mass function, such that the 
predicted number of clusters in the i'th bin would be 



where Vbcs ( ^m , ^^A = 1 — ^m) is the BCS survey volume 
defined in § 13.11 For the theoretical mass functions, we 
calculated the predicted number of observed clusters in 
four mass bins (of widths 5, 5, 5, and 10 x IQ-^^H^^Mq, 
respectively), ranging from 2.3 x lO^^/i^^A^© < Afigoc < 
7.3 X 10"/i-iMq to 17.3 x IO^/i-^Mq < Migoc < 
27.3 X IO^'^H'^Mq. In the end, we only used the three 
uppermost mass bins, since the uncertainties of the low- 
est mass bin will be very large, and completely domi- 
nated by uncertainties in the incompleteness correction, 
given the small degree of sample completeness in this 
mass range. In the lowest bin, the incompleteness cor- 
rection will be highly sensitive to the scatter around the 
mean Migoc — Lx relation. As shown in § 13.31 the sam- 
ple selection cutoff in Lx corresponds to a lower cutoff 
in mass of Migoc ~ 5 x 10^'^h~^MQ. To account for 
Poisson shot noise, ^/N errors are applied to the counts 
of clusters in each mass bin. Additional uncertainties 
arise from the cluster mass measurement uncertainties, 
resulting in a non-negligible probability that a cluster 
is counted in the wrong bin. These errors are added in 
quadrature to the Poisson errors. In addition, ^/N errors 
are applied to the theoretically predicted cluster counts 
in each bin, to account for cosmic variance. Fitting to 
the theoretical mass function of Sheth & Tormen (1999), 
we find cr8(rJm/0.3)0-37 = O.Qltom (68% confidence lim- 
its) , while a similar fit to the mass function of Warren et 
al. (2006) gives agi^m/O.Sf-^'^ = 0.66 ± 0.04. 

There is a tendency for a preference for higher erg values 
in the higher mass bins; the uppermost mass bin alone 
gives a best fit dg = 0.72toio for = 0.3, using the 
mass function of Sheth & Tormen (1999). This may hint 
towards an overestimate of the sample completeness in 
the lower mass bins. Stanck et al. (2006) discusses the ef- 
fect of Malmquist bias on X-ray flux-limited cluster sam- 
ples in the presence of large intrinsic scatter around the 
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AI — Lx relation and suggest that the true normahzation 
of this relation is a factor of 2 dimmer than estimated 
from the RASS-based HIFLUGCS sample of Reiprich & 
Bohringer (2002). Applying a shift of this magnitude to 
the Misoc — Lx relation derived in § 13.31 (while keeping 
the estimated level of scatter, which is consistent with 
the scatter estimated by Stanek et al. 2006) would in- 
crease the value of M(Lmin) (to 8.7 x lO^'^h'^MQ, in the 
case of the concordance model) and increase the incom- 
pleteness correction factors in each mass bin, yielding a 
best fit 0-8 (f^m = 0.3) = 0.74;o;og, based on the three 
uppermost mass bins. However, since the sample of 50 
clusters used to derive the Afigoc— -^x relation in S I3.3l can 
be more accurately described as a volume-limited rather 
than flux-limited sample, this Malmquist bias correction 
should not be applicable to our sample. 

It should be noted that the magnitude of the uncer- 
tainties in the gravitational lensing mass measurements 
of the clusters, including the three additional sources of 
mass scatter discussed in § 13.21 has a significant effect 
on the derived cosmological parameters (see e.g., Lima 
& Hu 2005 for a discussion of a similar effect on dark 
energy constraints from high-redshift clusters). The ef- 
fect of the scatter is to populate the high-mass tail of 
the cluster mass function with some of the (far more nu- 
merous) lower-mass clusters which have been scattered 
upwards in mass. As an example, if we artificially reduce 
the scatter in mass in our data by 30% (as would be the 
effect of ignoring the effect of mass scatter from varia- 
tions in the NFW concentration parameter, variations in 
cluster galaxy contamination and projection effects), the 
best fit value of erg for the concordance model would in- 
crease from as = O.GT^q'^I to as — Hence, a 
decrease in the mass uncertainties works in the opposite 
direction as an decrease in the mass values, so that these 
will partially cancel each other out if the cluster masses 
and mass uncertainties are multiplied by the same factor. 
This makes our erg result fairly robust to e.g. a modest 
multiplicative error in the shear estimator. 

4. DISCUSSION 

The results presented in this paper demonstrate that 
weak gravitational lensing based measurements of cluster 
abundances provide competitive constraints in the as-^m 
plane. In a recent work (Pedersen & Dahle 2006), as was 
estimated by measuring the normalization of the mass- 
X ray temperature relation of galaxy clusters. If the 
intrinsic scatter about the mass-temperature relation is 
low (10% or less), as this normalization implies a higher 
value of as = 0.88 ± 0.09 for flm = 0.3. While most the- 
oretical studies predict a rather tight mass-tempcrature 
relation, Pedersen & Dahle (2006) showed that the scat- 
ter about the relation is most likely larger than 10% , 
producing a cg estimate which is systematically biased 
towards high values. Hence, as remarked in their discus- 
sion, the CTg estimate of Pedersen & Dahle (2006) should 
most accurately be regarded as an upper limit, and the 
true value is likely to be lower. 

Being based directly on weak gravitational lensing 
mass estimates, the results presented here do not de- 
pend on knowing the mass-temperature relation and its 
scatter, and we obtain a lower value of cg — 0.67j^Q Q5 
for fljn = 0.3. As shown in Figure |31 our joint constraints 
on CTg and ilm agree very well with the 3-year WMAP 



data (Spergel et al. 2006). The 3-year WMAP data give 
a preference for as and f2m values lower than indicated 
by the first year WMAP data, combined with large scale 
structure data from SDSS (Tegmark et al. 2004) . Our re- 
sult also agrees well with many estimates of cluster abun- 
dances based on X-ray data (e.g., Reiprich & Bohringer 
2002; Seljak 2002; Viana, Nichol, & Liddle 2002; Allen 
et al. 2003; Henry 2004), which have tended to produce 
CTg values in the lower end of the range of published es- 
timates. Reiprich (2006) discusses cluster abundances 
in light of the 3-year WMAP results, and gives best-fit 
results for CTg and 17™ based on the HIFLUGCS cluster 
sample that agree very well with both WMAP data and 
the results presented here. These results imply a lower 
value of CTg than implied by "cosmic shear" measurements 
of weak gravitational lensing in random fields (Hoekstra 
et al. 2005; Semboloni et al. 2006; see figure|21). 

One possible source of systematic error in the results 
presented in this paper is the use of X-ray luminosity 
selected clusters to estimate the scatter about the mass- 
X-ray luminosity relation. Some weak lensing studies 
(e.g., Erben et al. 2001; Dahle et al. 2003) may hint to- 
ward the existence of a class of X-ray under- luminous 
clusters, although detections of shear-selected clusters 
are significantly affected by projection effects (White, 
van Waerbeke, & Mackey 2002; Hennawi & Spergel 
2005), and may include "proto-clusters" (Weinberg & 
Kamionkowski 2002) that are still undergoing collapse 
and have not yet become fully virialised systems. The 
latter type of objects may pose less of a problem for the 
present study, since such systems should not be included 
when making the fit to a theoretically predicted mass 
function for fully collapsed systems. A sub-population of 
X-ray under- luminous clusters is also found by Popesso 
et al. (2006), based on a correlation of Abell clusters 
in the SDSS database and RASS sources. Since these 
systems show signs of being in formation, still undergo- 
ing significant mass accretion, their relation to theoreti- 
cally predicted abundances of virialised clusters is some- 
what unclear. Weak lensing surveys with the large sky 
coverage required to find statistically useful numbers of 
Misoc > IO^^Mq clusters are still some years into the 
future, but ongoing surveys (Miyazaki et al. 2003; Het- 
terscheidt et al. 2005; Gavazzi & Soucail 2006; Wittman 
et al. 2006; Schirmer et al. 2006) are beginning to de- 
tect a significant number of lower-mass clusters based on 
weak gravitational shear. If there is indeed a large popu- 
lation of massive, X-ray under-luminous galaxy clusters, 
our estimates of sample incompleteness as a function of 
mass would be too low, implying a somewhat higher true 
value of CTg. 

The constraints presented here may be improved in 
several ways in the near future: Firstly, obtaining new 
wide-field imaging data for all the clusters in this sam- 
ple will remove the need for extrapolation of an adopted 
mass model out to ^igoc, removing one potential source 
of systematic error while significantly reducing the statis- 
tical errors on the mass measurements. Furthermore, by 
obtaining imaging in a larger number of pass-bands, ap- 
proximate photometric redshifts may be calculated, pro- 
viding sufficient 3D information to distinguish the grav- 
itationally lensed background galaxies from cluster and 
foreground galaxies. Although we have made an attempt 
to correct for this effect by applying an average, radi- 
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Fig. 3. — Joint constraints on trg and Qm- The blue curves show 68% (thick Unes) and 95% (thin lines) constraints on these two 
parameters from the galaxy cluster data presented in this paper. Also shown are constraints on the same parameters from the CFHTLS 
weak lensing survey, from the 3-year WMAP CMB data, and joint constraints from these two data sets, derived by Spergel et al. (2006). 



ally dependent, contamination correction to all clusters 
(and estimate the extra scatter caused by cluster rich- 
ness variations), the true correction factor will of course 
depend on cluster richness and cluster redshift, and clus- 
ter substructure will cause it not to be a simple function 
of cluster radius. Hence, the appropriate correction for 
individual clusters can differ significantly from the av- 
erage correction calculated here (see e.g., Broadhurst et 
al. 2005). Future photometric redshift estimates of back- 
ground galaxies, and the uncertainties in these estimates 
could also be included in a optimized weighting scheme 
(e.g.. Kaiser 2000; Benitez 2000), significantly improv- 
ing the signal-to-noise ratio of the lensing mass measure- 
ments. 

The sample used in this study is restricted to galaxy 
clusters in the Northern celestial hemisphere, and the 
sample size could easily be doubled by including clusters 
with similar X-ray selection criteria in in the Southern 
hemisphere, e.g., from the REFLEX sample of Bohringer 
et al. (2004), thus improving the statistical errors. The 
sample could also be extended by lowering the X-ray lu- 
minosity cutoff Lyaim Sampling the mass function at lower 
masses. For the (e)BCS sample, such a low-luminosity 
extension of the cluster sample would be restricted to 
the near end of the redshift range probed here, since 
■^min corresponds to the lower eBCS X-ray fiux limit at 
the far end of the survey volume. However, more re- 
cent cluster surveys based on RASS data (e.g., Ebeling, 
Edge, & Henry 2001) go down to lower X-ray fluxes and 
would be useful for probing further down the mass func- 



tion without restricting the survey volume. Finally, by 
combining the data presented in this paper with a weak 
lensing survey of a similar cluster sample at higher red- 
shifts, interesting constraints can also be put on the dark 
energy equation of state parameter w, based on the evo- 
lution of the cluster mass function. 
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